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O _ 

O ; Abstract 



We present several second-order linear differential equations that are associated to a particular 
Riccati equation with only one constant parameter in its coefficients through the technique of 
Q ! supersymmetric factorizations and through a Dirac-like procedure. The latter approach is a 
^ I minimal extension of the results obtained with the first technique in the sense that it includes 
^ ■ up to two more constant parameters. 

^ ; PACS Numbers: ll.30.Pb 

o 

^ ! 1. Riccati solutions and factorizations 

(N . 

O ■ We start with one of the simplest Riccati equation, which nevertheless dictates for example the 
^ ' behaviour of the Hubble constant of barotropic Friedmann-Robertson-Walker cosmologies [ 1 ] 

' where fc = ±1 and c is a real constant. Employing u = one gets the well-known second 
order differential equation 

wl+Kc^w^ = 0. (2) 

' For K" = 1 the solution of the latter is Wj = Wj cos(c77 + (p), where (p is an arbitrary phase, 
whereas for k = —I one gets w j = W jsinh(c77) where W_j_j are amplitude parameters. 

The point now is that the Riccati solution Up = ^ ^ mentioned above is only the particular 
solution, i.e., Up = — tan(c77) and Up = coth(cT]) for K = ±1, respectively. The particular 
Riccati solutions enter as nonoperatorial part in the common factorizations of the second-order 
linear differential equations that are directly related to the well-known Darboux isospectral 
transformations [ ]. Indeed, Eq. (2) can be written 



u' + cu^-FK-c = , (1) 



w -c(-K-c)w = (3) 
and also in factorized form using Eq. (1) one gets (D^ = ^) 

(D^ + cUp) (D^ - cUp) w = w" - c(Up + cu2)w = . (4) 

To fix the ideas, we shall use the terminology of Witten's supersymmetric quantum mechanics 
and call Eq. (4) the bosonic linear equation. On the other hand, the supersymmetric partner 
(or fermionic) equation of Eq. (4) will be 

(D^ - CUp) (Dj, + CUp) Wj = Wf - c(-Up + cup)wf = Wf - c • c^ ^ (77) w^ = , (5) 
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which is related to the fermionic Riccati equation 

-u'+cu2-c^f(r7) =0 . (6) 

Notice that for this fermionic Riccati equation the free term is already not constant. Explicitly, 
one gets 

c(l +2tan2c77) if K" = 1 
c(-l +2coth^c77) ifK- = -l 



-Up + CUp 



w 



metric j 

for K = 1 and K" = — 1 , respectively. 



for the supersymmetric partner Riccati free term. The solutions are w^^ = ^^^^^'^^^■^ and 



f sinh(cTj 

Introducing the operator = — zDj^ we can write the fermionic equations as follows 

(-P^ -icup) (Pjj -icup) Wf = -P^w^-c(-iP^Up + cUp)Wf 

whereas the bosonic case is 

(P^ -icUp) (-Pjj -icUp) Wj, = -p2wj,-c(iP^Up + cu2)w;, , 

There is a more general bosonic factorization that has been introduced for the case of the 
quantum harmonic oscillator by Mielnik [5] that in our case reads 



(7) 



(8) 



Djj +cUgfj [Dri -cUgf j Wg = Wg -c(Ug f + cugf)wg = Wg + K-cc(77; A)wg = .(9) 

It is given in terms of the general Riccati solution Ug ^(t]) of the fermionic Riccati equation 
(6) 



Ugf(T7;A) =Up(77)--D^ ln(I^(77)+A) 



In 



(10) 



and yields the one-parameter family of Riccati free terms c^{ri;X] 



-Kc^(77;A) =cuJ(77;A) + 



dug(r];A; 



drj 



ln(I^(7])+A: 



jv^MKil). 2w4(^) 



(11) 



c(Ut7) + A) c(I,(t])+A)2' 

where Ik{v) = Iq w^(y)dy, if we think of a half line problem for which A is a positive 
integration constant that occurs as a free parameter of the method. 

The free terms 0^(7]; X) have the same supersymmetric partner free term c^^{r}). They 
may be considered as intermediates between the initial constant Riccati free term kc and its 
supersymmetric partner free term c^^{r}). From Eq. (9) one can infer the new parametric 
'zero mode' solutions of the linear equations entailing the functions c^{r};X) as follows 

For A ^ oo the solutions Wg vanish. For an application to the damped classical oscillator 
see In quantum mechanics, where we have different types of solutions and physical 
interpretation, one can obtain Wg in the same asymptotic limit of the parameter if one 

introduces the appropriate normalization constant [-f] . 
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2. Dirac-like approach 

The Dirac equation in the susy nonrelativistic formalism has been discussed by Cooper et al 
[6] already in 1988. They showed that the Dirac equation with a Lorentz scalar potential is 
associated with a susy pair of Schroedinger Hamiltonians. This result has been used later 
by many authors in the particle physics context. In mathematical terms, it is only a simple 
approach for matrix differential equations. Here we present several applications that may be 
considered the simplest extension of the results of the previous section in the sense that only 
up to two more constant parameters are introduced. 

2.1 - Let us now consider the following two Pauli matrices a = —iOy = ~i ^ ^ 
and /3 = Of = ^ ^ q ^ and write a Dirac-like equation for zero mass and at fixed zero energy 



of the form 



Dl = [iayP^ + ax(icup)]W = , (13) 



where ^ ~ J ^ '^^^ component 'zero-mass' spinor. This is equivalent to the 

following decoupled equations 

— Pj^Wj -|-icupWj = (14) 
+Pj^W2 + icUpW2 = . (15) 

Solving these equations one gets Wj oc 1/ cos(cT]) and W2 °^ cos(cr7) for the K = \ case and 
Wj l/sinh(c77) and W2 °= sinh(c77) for the K" = — 1 case. Thus, we obtain 

w=(Z'] = (:/]. (16) 



This shows that the matrix 'zero-mass' Dirac equation is equivalent to the two second-order 
linear differential equations of bosonic and fermionic type, Eq. (2) and Eq. (5), respectively 
[6]. 

2.2 - Consider now the following Dirac equation 

D2 = [iayP^ + ax(icup + K)]W = KW , (17) 

where ^ is a positive real constant. In the left hand side, K stands as a mass parameter of 
the Dirac spinor, whereas on the right hand side it corresponds to the energy parameter, i.e., 
E = K. Thus, we have a Dirac equation for a spinor of mass K at the fixed energy E = K. 
This equation can be written as the following system of coupled equations 

-P;,Wi-^(icUp-FK)wi =Kw2 (18) 
PT]W2 + (icUp-FK)w2 = Kwi . (19) 

This system is equivalent to the following second order equations for the two spinor 
components, respectively 

i(TP^-2K)up + cu2jw^ = 0, (20) 



where the subindex i = 1,2. 

The fermionic spinor components can be found directly as solutions of 



D2 W+ 



c^(H-2tan2c7]) + 2icKtanc77 W]^ = forK-=l (21) 
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and 



c^(-l +2coth^cr]) -2icKcothcT] 



w, 







whereas the bosonic components are solutions of 



and 



D>2 + 



D^W2 + 



c — 2icKtanc77 







•c +2icKcothc77 







for K=l 



for K 



forJC : 



1 . 



(22) 



(23) 



(24) 



The solutions of the bosonic equations are expressed in terms of the Gauss hypergeometric 
functions in the variables y = e''^'' and y = e'^^, respectively 



_ 1 i 2 



1 

2' 



+ 



2F1 (p - iq) , ^ (p + iq) , 1 + p; -y^ 



(25) 



and 



W2 =C(-l)-^V"2Fi[-^(r + is),-^(r-is),l-ir;y2 



+ 



D( 



-i)^y%F 



2^ 1 



^(r-s),^(r + s),l+ir;y2 



(26) 



The parameters are the following: p = ( — 1 — q=(l 
3, s = ( — 1 +i— whereas A, B, C, D are superposition constants. 



respectively. 

r=(_i_if ; 

It is not necessary to try to find the general fermionic solutions through the analysis 
of their differential equations (21) and (22) because they are related in a known way to the 
bosonic solutions [ /]. The general fermionic solutions can be obtained easily if one argues 
that the particular fermionic zero mode is the inverse of a particular bosonic zero mode and 
constructing the other independent zero mode solution as in textbooks. Thus 

^_l+kf[w2^fdz 



W, 



w. 



± 



(27) 



where k is an arbitrary constant. 



2.3 - The most general case in this scheme with only constant parameters is to consider 
the following matrix Dirac-like equation 



D3 = 



-i 

1 



Pt7 + 



1 

1 



icUp + Kj 
icug + 



Wj 
Wn 



Kj 




Wj 

VJn 



Proceeding as in 2.2 one finds the coupled system of first-order differential equations 



+ icug + K2 



Wn 



KjWj 



Pjj +icUp + Kj 



w, 



K2W2 



(28) 

(29) 
(30) 
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and the equivalent second-order differential equations 



ic(up-Ug) + (Kj-K2) 



ic(±P7^u. + K^Ug + K2Up) -c^UpUg 



w.. = 



(31) 



where the subindex i = 1,2, and u^ and Uj correspond to Up and Ug, respectively. In the D,^ 
notation this equation reads 



cAupg - iAK) 



+ 



c(±Dj^u.-|- (iK^Ug-|-K2Up)) -c^UpUg 



w.. = . 



Under the gauge transformation 

1 



w,. = z.exp 



cAupg-iAK) dz] =z,(t]) 



e3 



iiTjAK 



(I^ + A)2 



one gets 



where 



-D„ 
2 



Prih^ + Q,{ri)z, = 0, or D^\. + 2,(77)2, = 0, 

QiiV) = c(±D^u,. + (iKiUg + K2Up))-c-UpUg 
cAu, 



(32) 



(33) 



(34) 



1 r i2 

- cAupg-iAK) 



(35) 



for / = 1,2, respectively. Since 2- are complicated functions we were not able to find 
analytical solutions of Eq. (34). 

The corresponding Dirac spinor is of the following form 



W(7];Ki,K2,A: 



2i(t?;Ki) 

z2(^7;K2,a: 



Wf(T];Ki) 
Wg(77;K2,A: 



(36) 



where Wg(77;K2, A) is givenby Eq. (12) whenKj = K2 = 0. The discussion of the asymptotic 
limit A ^ 00 is similar to that at the end of the first section. 
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